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Abstract 

If the neutrinos are to be identified with the primary source of ultra-high energy cosmic 
rays(UHECR), their interaction on relic neutrinos is of great importance in understanding 
their long intergalactic journey. In theories with large compact dimensions, the exchange of 
a tower of massive spin-2 gravitons (Kaluza-Klein excitations) gives extra contribution to 
vv — > // and 77 processes along with the opening of a new channel for the neutrinos to 
annihilate with the relic cosmic neutrino background vv — > Gkk to produce bulk gravitons 
in the extra dimensions. This will affect their attenuation. We compute the contribution of 
these Kaluza-Klein excitations to the above processes and find that for parameters of the 
theory constrained by supernova cooling, the contribution does indeed become the dominant 
contribution above -Js — 300 GeV. 
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Ultra high energy (UHE) cosmic rays [[3J present a challange in cosmology and 
astroparticle physics. Energies beyond the Greisen-Zatsepin-Kuzmin (GZK) cut 
off ~ 5 x 10 19 eV excludes protons as primary candidates || for high energy 
events due to the fact that above the resonance threshold for A* production, the 
protons would lose roughly 20 % of their energy by scattering on the 2.73 K 
cosmic background photons through p + j 2 . 73K — ► A* — > N + % during 6 Mpc 
(mean free path) of their journey. Thus apart from the difficulty in accelerating 
particles to such extremely high energies 0, these protons would lose significant 
fraction of their energy if they came from 100 Mpc or further away. Gamma rays 
and neutrinos could be other primary candidates which will not be hampered by 
the GZK cut off. The high energy photon annihilation cross-section on cosmic 
background photons Q will however, preclude them from being the candidates 
particularly if they have to travel a distance of the order of 100 Mpc of the 
intergalactic space, their mean free path being of the order of 10-40 Mpc only. 
Weiler [|[ proposed that if neutrinos have to be the source of UHE cosmic rays, 
they have to have their energies corresponding to the Z-resonance in order to be 
able to efficiently annihilate relic neutrinos i.e. E v ~ |^ ~ 10 23 eV for neutrino 
mass consistent with the Super-Kamiokanda data which is clearly above the GZK 
cut off |J. The difficulty in realizing this scenario is to identify the source of ultra 
high energy neutrinos with their energy close to the Z-resonance. Decay of a super 
heavy relic particle into neutrinos with mass M x > 10 13 GeV has been proposed 
H to be the source of these highly energetic neutrinos that can explain cosmic ray 
events beyond the GZK cut off. Such neutrinos however, have to be produced in 

a spherical shell at the red shift z = {j^r- 1) = { M ^^ — 1) so that their energy 

near earth is close to the Z-resonance energy for the vv annihilation cross-section 
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on the relic neutrinos to be large. If such fine tuning of mass of the relic super 

heavy particle or their abundance at the red shift z = {^e^ 1) is not available, 

we will require large neutrino hadron cross-sections of hadronic strength for the 
neutrinos to be able to initiate showers high in the atmosphere. Current estimates 
of the ultra high neutrino scattering in the Standard Model for E v ~ 10 21 eV 
are estimated [0] to lie in the range 10~ 4 — 1CT 5 mb. Specifically, 

a uN SM (E u ) * (2 - 3)xlO- 5 ( T J^ F )°- 4 mb (1) 

The anomalously large high energy neutrino interaction requires new physics 
and an appealing possibility considered in the literature is the theory with n extra 
dimensions with large compactification radii and TeV scale quantum gravity ||. 
In such theories, exchange of a tower of massive spin-2 bulk gravitons (Kaluza- 
Klein(KK) excitations) opens up the exciting possibility of extra contribution to 
any two particle scattering. In this context, the contribution of gravitational 
scattering to vN cross-section consistent with the requirement of unitarity has 
been estimated @. For example, for a neutrino of energy E v scattering on a 
proton at rest 

4-7TS \TeV E 

^) = ^<W^ (2) 

The most stringent lower limits on the string scale Ms are obtained from astro- 
physical considerations |10| by studying the energy loss rates through the emission 
of bulk gravitons in SN 1987 A and lead to typically Ms > 30 Tev, 4 TeV and 1 
TeV for the number of extra dimensions equal to 2, 3 and 4 respectively. Thus 
n = 4 is required if neutrinos are to be the primary candidates for UHE cosmic 
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ray events. 

The new physics would open up new channels for these high energy neutrinos to 
annihilate with the relic neutrino background to produce bulk gravitons in the 
extra dimensions through 

vv — ► G KK (3) 

and may therefore get effeciently attenuated during their intergalactic journey. 
In these theories, there also exists extra contributions to high energy neutrino 
scattering on relic background neutrinos through the exchange of KK excitations 
to produce the cascade of high energy particles through 



— > // (4) 

and 



vv 



77 (5) 



processes. In the Standard Model the process (4) takes place through the exchange 
of Z, while the process (5) can take place only through the loop diagrams and is 
expected to be highly suppressed Jllfl . In theories with extra dimensions, there 
is direct coupling between Standard Model particles and spin-2 gravitons, thus 
providing a unique channel. 



Using the Feynman Rules given in |TJ] , it is straight forward to calculate the above 
processes. The spin averaged neutrino annihilation cross-section to produce a KK 
graviton can be easily found to be 

7T K/^ f~S 

oivv — > G K K{s,m n )) = — — — S (m n - y/s) (6) 



4 



where m n is the mass of the nth KK state, being given by m} n = 47 ^" 2 where 
k 2 R n = 16tt (4tt)^ r(| ) M~ {n+2 \ Summing over all the KK states using the 



density of the number of states given in equations (B2) and (B3) in []12 , the total 
cross-section comes out to be 

a{ V u^G KK ) = £(^)*« (7) 
where n is the number of extra dimensions. 

The matrix elements for the scattering processes (4) and (5) are given in |T2j (see 
equations (61), (62), (71) and (72)). Squaring the matrix elements and summing 



over the virtual KK state exchange propagator using equations (B4)-(B8) of JT2 
the total scattering cross-section can be calculated in a straight forward manner 
and we get 



_ 7T S n+2 

oHvv — > ff) = (-^) T 2 (8) 



and 



7T « n+2 

oHvv — > 77) = (Ay) T 2 (9) 

where T 2 = tc 2 + 4/ 2 (^j) and / is given in |T2| (see equation (B8)). 
In these calculations it was assumed that the massless four dimensional graviton 
and its massive KK excitations couple with the same strength and the sum over 
all KK states to a given scattering amplitude in general is obtained by introduc- 
ing a cut-off. However, it has recently been pointed out that due to brane 



fluctuations, effective coupling of the nth KK mode to four dimensional fields is 
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suppressed exponentially and is given as 



— c vnt 



K n = KV^\—^\ (10) 



s 



where c parametrises the effects of finite brane tension and is of order one or 
larger. This suppression then provides a dynamical cut-off in the sum over KK 
states. With this modified coupling, neutrino-nucleon cross-section due to the 



exchange of KK gravitons has been recalculated [H and has been found to be 
considerably smaller to account for UHE cosmic ray observations. If that is so, 
neutrinos remain deeply penetrating and can be the source of UHE cosmic rays 
only if their energies correspond to the Z resonance 0,0. 
For a neutrino of energy E u annihilating a relic neutrino of mass m v , we get 



„ s 3 „ »/m„\a/ K \ 2 fl'l'eV\ n+2 9 , , 

°» - G«rf - 4 x ID"— 2, (^) (^) (_) cm' (11) 

and 

(12) 

and 

o 9 {yv — ► 77) = 3a 9 (uu — > ff) (13) 
to be compared with the Standard Model calculation 

a su {v - _^ fJ) = G% Mis 2 2 

where C v f = T 3 f sin 2 9 w and Cy = . The Standard Model cross-section 
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at s — Mj comes out to be 



1.59x10 cm . The process vv — > 77 has 



been recently calculated in the Standard Model [15] and is given by 



a SM (vv 



77) 



G 2 F a 2 s 3 
640^ 3 



A 2 ~ 5.6 x 10" 42 



m,. 



leV 



10 2 W 



cm 



(15) 



where A is taken to be ~14.4. The above expression is valid for yfs upto roughly 
2Mw after which the cross-section reaches a plateau till about 1 TeV and then 



starts falling [15 



The effect of brane fluctuations through modification of the coupling can be easily 
incorporated in the calculations. a{yv — > GkK) gets multiplied by a suppression 
factor exp( z ^) and in the scattering cross-sections (8) and (9), T 2 is replaced by 
\T>\ 2 where 



V 



— ltv exp 



2cs 
M 2 



- J 



2cs 
M 2 



with 



J 



2cs 



V \ dy- exp 

Jo 1 - y 



-2cy 



(16) 



(17) 



Ml) Jo a l - y ^\ Ml J 
where V denotes the Principal Value of the integral and for even values of n it 
reduces to 



J 



2cs 

Ml, 



1 



,71 — 1 



Ei 



2cs 
Ml, 



2 — 1 
2 L 



+ Y, (k- 1 ) 1 ex P 

fc=i 



2cs 
Ml, 



Mi 



JT 

2cs 



18) 



where Ei{x) is the Exponential-integral function [16| . 

The suppression in the cross-sections due to brane fluctuations at the available 
values of yfs even for the smallest allowed value of M s = 1 TeV for n = 4 is not 
more than 10 % and is not significant as can be seen from the expression of T>. 
In Figures 1-3, we have plotted the cross-sections for the cases n = 2 and 4 with 
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M s = 1, 10 and 30 TeV. In Fig.4, we plot the cross-sectons for the parameters of 
the theory that are constrained from astrophysical considerations obtained by de- 
manding that the supernova cooling through graviton production does not exceed 
10 53 ergs. We have also plotted the Standard Model cross-sections for vv — > // 
and vv — > 77 for comparison. We see that the contribution of spin-2 graviton 
exchange for n=2 and M s = 1 TeV rises very fast with energy and far exceeds 
the Standard Model contribution. Such a large cross-section would result in the 
attenuation of high energy neutrinos if they are produced through the decayof 
a super heavy relic particle of mass M x exceeding 10 14 GeV. The probability of 
these high energy neutrinos to interact depends upon the cross-section and on 
the lepton-asymmetry parameter 77 = n "^ el ™ whose value is constrained by nu- 
cleosynthesis and by large scale structure considerations. The mean free path of 
neutrinos is dominated by oivv — > Gkk) and is given by A = —t-t. — — « . 

U\UV ■ KK )Ylv, relic 

If there is no lepton-asymmetry, 77 ~ 0.14 and could be as large as 4 for some 



cosmological parameters fL7 |. For y/s = 500 GeV,the neutrino mean free path 
for n = 2 and M s = 1 TeV is \ u v^g kk ^ 0.57xl0 7 (^) Mpc to be compared 
with the corresponding contribution in the Standard Model from the Z-exchange 
namely \ u9 ^z* — 1.33x 10 10 (^ 1 ^) Mpc. For the astrophysically prefered value 
(n = 4, M s = 1 TeV), the mean free path is X u d-,g kk - 2.27 x 10 7 (^) Mpc. 
At resonance, y/s = Mz, the Z-exchange gives by far the largest cross-section with 
mean free path \ vi ^z* Te , a — 0.345xl0 6 (^^) Mpc to be compared with 
\ u &->Gkk — 6.8xl0 8 (^^) Mpc for the astrophysically allowed parameters. 
However, neutrinos produced by the decay of super heavy relic particles with 
yfs > 500 GeV would mainly annihilate to produce bulk gravitons with a cross- 
section more than three orders of magnitude greater then the Standard Model 



s 



cross-section through Z* exchange. It should however be noted that these new 
neutrino interactions with CMB photons and massless relic neutrinos would not 
have enough center of mass energy to compete with the Standard Model interac- 
tions to significantly influence the neutrino propagation. For massive neutrinos 
on the other hand, the relic neutrino background provides large center of mass 
energy for these interactions to play an important role in neutrino propagation 
through their intergalactic journey [18|]. For example, from Fig. 4, we observe 
that for Super-Kamiokande motivated neutrino mass ~ 0.07 eV, vv annihila- 
tion cross-section for bulk graviton production for neutrinos of energy ~ 10 15 
GeV, (y/s ~ 265 GeV) is comparable to the Standard Model cross-section. It 
is thus clear that in a scenario in which the neutrinos are produced through the 
decay of super-heavy relic particles of mass Mx > 10 13 GeV, the vv annihilation 
cross-section on the relic neutrinos in the Standard Model is large enough at ^fs 
near the Z-resonance energy to qualify for the source of UHE cosmic rays beyond 
the GZK cut-off even for Super-Kamiokande constrained neutrino mass range. At 
these energies the new interactions remain much smaller than the Standard Model 
ones. This of course, requires that the X particles have life time long compared 
to the age of the universe, Tx > 10 5 t Q . On the other hand, if neutrinos were 
produced at large z i.e. at large cosmological distances, which would be the case 
for super-heavy particles with life time much smaller than the age of the universe, 
large center of mass energy (enhanced by a factor of (1+z)) will be available to the 
neutrinos for new interactions to play an important and for values of y^i > 265 
GeV even dominant role. The cascade of charged particles and 7 rays produced 
through the exchange of a tower of massive KK excitations in theories with n ex- 
tra dimensions with large compactification scale will give significant contribution 
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to the GeV gamma ray background today. 
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where in they have proposed that from the measurements of horizontal air 
showers one can put bounds on the UHE neutrino flux striking the earth 
which in turn can be used for putting bounds on UHE neutrino interactions 
with matter beyond the Standard Model. Thus the modification in neutrino 
propagation in theories with large extra dimensions can in principle be tested. 
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Total vv annihilation cross-section as a function of y/s. Curves a, b, c, d and f stand 

G KK ,n = 2), o°{vv — > ff,n = 4), 



M (yv — > //), — > //, n - 2), - 

* GxK,n = 4) respectively for M s = 1 TeV. 
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Figure 3: Total vv annihilation cross-section as a function of y/s. The curves are as in Fig.l 
and M s = 30 TeV. 
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Figure 4: Total vv annihilation cross-section as a function of yfs for the astrophysically 
trained values of the parameters namely, M s — 30 TcV for n = 2 and M s — 1 TeV for n 
Curves a, b, c, d and f are as in Fig.l. Curve g stands for o SM {vv — ► 77) 
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